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We show how hot QCD equations of states can be adapted to make definite predictions for 
quark-gluon plasma at RHIC. We consider equations of state up to 0{g^) and 0[g^ {ln{l/ g) + 
5)]. Our method involves the extraction of equilibrium distribution functions for gluons and 
quarks from these equations of state by capturing the the interaction effects entirely in the 
effective chemical potentials. We further utilize these distribution functions to study the 
screening length in hot QCD and dissociation phenomenon of heavy quarkonia states by 
combining this understanding with the semi-classical transport theory. 
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I. INTRODUCTION 

The hot and dense matter created at RHIC behaves like a near perfect fluid as inferred from 
recent experimental observations at RHIC[l[. This form of the matter lies in the strongly inter- 
acting domain of QCD. It is of the interest to check if the domain is genuinely non-perturbative, 
or could be understood in the form of higher order terms in perturbation theory. Methods based 
on weak perturbative techniques 0] have achieved partial success in interpreting the results from 
RHIC. The aim of the present article is to investigate more fully, the viability of pQCD EOS in 
understanding the RHIC results. We implement this programme by extracting quasi-free equi- 
librium distribution functions from the EOS , capturing all the interaction effects in the effective 
fugacities[^, 13]. This procedure allows one to study various observables for hot and dense matter 
at RHIC. 



II. QUASI-PARTICLE REALIZATION OF HOT QCD EOS 



We consider hot QCD EOS in pQCD up to 0(g5)[3](EOSl) and 0[g^{\n[l/ g)+5)\^{F.0^2). The 
free parameter 5 in E0S2 encodes the unknown perturbative and non-perturbative contributions 
of 0{g^) and can be fixed by matching the EOS with lattice predictions on EOS[H|. The ansatz 
for the determination of equilibrium distribution functions involves retaining the ideal distribution 
forms, with the effective fugacities {zg = exp{iig) and Zf = exp(/x/)). Note that for the mass less 
quarks ( u and d) which we consider to constitute the bulk of the plasma, Zg/ / = 1 {fJ-g/ / = 0) if 
they were not interacting. 

As the first step in our approach, we express the hot EOS in the form 

P = P^+P^ + APg+APf (1) 

We seek to absorb the contributions from all the non-ideal terms in the effective chemical potentials 
/ig and fif. Since the equations of states are expected to be valid at T > 2Tc, we treat the 
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dimensionless quantity jlgj = PfJ.gj perturbatively. This determination of p-g/f needs to be done 
self consistently. Accordingly, we expand the the grand canonical partition functions (Zg^f) for 
gluons and quarks as a Taylor series in p-gj. We obtain 



fe=0 



fc=0 



(2) 



We equate log(2') with the pressure by the well known relation : log(Z) = P/3V, and determine 
pg/f self-consistently order by order. We find that it is sufficient to determine them up to cubic 
order. Details may be found in 



III. THE SCREENING LENGTH 

The Debye mass from equilibrium distribution function can be obtained from the expression 



Mt 



■'g/f 



i9')2Cg/f 



d f^/f 



(3) 



where Cg = 2Nc, Cf — —"^Nf and g' is the effective coupling which appears in the transport theory. 
We fix g' by comparing the temperature dependence of the screening length with the lattice results 
of Kazmarek and Zantowfioj. The behavior of the screening length for pure gauge theory and 
full QCD for EOSl and E0S2 are shown in Fig.l. Interestingly, the behavior of the screening 
length with temperature qualitatively matches with the lattice QCD results(see Fig. 2 of Ref . jTol| ) . 
The agreement for E0S2 is slightly better as compared to EOSl and in both the cases agreement 
becomes better at higher temperatures(T = 3Tc and higher). 




Pure DCD; S = 0.7 '■ 
Pme QCD; i = 0.9 - 




FIG. 1: Twice of the Debye screening length Ad is plotted as function of T/Tc for Pure QCD, Nf — 2and 
Nf — 3 for EOSl and E0S2. Please note that the phrase Pure QCD is used for pure gauge theory. The 
left Fig. is for EOSl and right one is for EOS2. Note that the screening length has been measured in 



IV. THE HEAVY QUARK POTENTIAL AND DISSOCIATION OF QUARKONIA 



In this section, we address the screening of heavy quark potential by combining the quasi-particle 
model for hot QCD introduced in Section(I) to the semi-classical transport theory @. The response 



function (chromo-electric permittivity) has recently been determined from the transport theory by 
employing the ideal distribution functions for gluons and quarks by Ravishankar and Ranjan[^. 
We determine these response functions for realistic distribution functions obtained from pQCD. 
Finally, the response functions is then used to study the screening properties of QGP. We outline 
the method(discussed in detail in[l, 0, @1) below. 

Consider the Cornell potential: (/)(r) = — -p+Ar, where a and A are phenomenological constants. 
The first term dominates at small distance while the linear causes confinement, dominating at 
large distances. The medium modifies the expression for the potential (in the Fourier space), 
^(fc) — > (f>{k)/e{k), where the response is evaluated at cj = 0. Accordingly, the potential undergoes 
modification 4){k) — > (j)m{k), which can be written as 

2 a 4 A 



where 



Nff2{Zq)+g2{Zg) 



(5) 



For a gluonic plasma, the above expression reduces to rn^jjg = l6nQ'^T^g2{zg), where Zg = exp(/Zg) 

and Zq = exp(/iq) Note that, mu has determined by the permittivity via: e = 1 + -jf-. The 
modified form of the potential in r-space is then given by: 

2A exp(-TOi5r) 2A 2A 

0s(r) = a) y 1 aniD- (6) 

mjj r mj^r mjj 

It is clear from Eq.(6) that the role played by the Debye mass in QGP is rather different from 
its role in electrodynamic plasma. To see this, note that at large T, Eq.(6) reduces to 

2A 

(l)s{r) ~ n aniD (7) 

Ignoring the additive contribution, the energy of the qq in the ground state is simply given by 

^9 - — (8) 

where rrig is the mass of heavy quark. The binding energy is, of course, temperature dependent 
and approaches zero as T — > cx). At any finite temperature though, the quarks possess a thermal 
energy Eth ~ f T, by equipartition theorem, leading to an ionization of the quarkonium when 
Etk matches the binding energy. This leads to the melting temperature Td of J/vP and T listed 
in Table (I). It is noteworthy that the dissociation temperatures are all roughly in the range Td « 
(2 — 3)Tc, which is higher than the temperatures achieved so far. Since the temperatures expected 
at LHC is ^ 2Tc — 3Tc, one may expect to test these predictions there. And moreover these 



estimates are consistent with the recent results from other theoretical works 11 



V. SUMMARY AND CONCLUSIONS 



In conclusion, we have developed a quasi particle model for hot QCD to extract the equilibrium 
distribution functions for gluons and quarks from the hot QCD EOS. We have shown that the 
interaction effects can entirely be captured in the effective fugacities for the gluons and quarks. 



TABLE I: The dissociation temperature(7D) for various quarkonia states (in unit of Tc). 



Hot EOS 


quarkonium 


Pure QCD 


Nf=2 


iV/=3 


EOSl 


J/* 


2.2 


2.62 


2.46 




T 


2.5 


3.14 


2.94 


EOS2 


J/* 


1.86 


2.38 


2.24 


S = 0.8 


T 


2.12 


2.76 


2.58 


E0S2 


J/* 


1.95 


2.45 


2.32 


5 = 1.0 


T 


2.2 


2.83 


2.66 


EOS2 


J/* 


2.03 


2.52 


2.40 


<5 = 1.2 


T 


2.28 


2.9 


2.74 



We utilized these distribution functions to study the screening length as a function of tempera- 
ture. We have determined the dissociation temperatures for J/vl/ and T by studying the medium 
modifications to heavy quark potential which enters in the form of chromo-electric permittivity. 
The results on dissociation temperatures are consistent with other recent theoretical works. Our 
approach can easily be generalized to the hot QCD EOS at finite quark-chemical potential. This 
approach can also be employed to study the thermodynamic and transport properties of hot and 
dense matter at RHIC. 
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